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Abstract 

In many models, stability of dark matter particles is protected by a conserved Z 2 
quantum number. However dark matter can be stabilized by other discrete symmetry 
groups, and examples of such models with custom-tailored field content have been pro- 
posed. Here we show that electroweak symmetry breaking models with TV Higgs doublets 
can readily accommodate scalar dark matter candidates stabilized by groups 7L P with any 
p < 2 N ~ 1 , leading to a variety of kinds of microscopic dynamics in the dark sector. We 
give examples in which semi-annihilation or multiple semi-annihilation processes are al- 
lowed or forbidden, which can be especially interesting in the case of asymmetric dark 
' ^ ' matter. 
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1 Introduction 

m 

Despite compelling astronomical evidence for existence of dark matter [TJ, there is still no 
direct experimental clue of which particle can play the role of dark matter candidate. It is 
only known that dark matter cannot be satisfactorily explained by the Standard Model (SM) 
particle content. In this situation, one can focus on exploring dark matter candidates arising 
in various models beyond the SM, especially if these models can simultaneously address other 
particle physics issues such as electroweak symmetry breaking and small neutrino masses. 

Dark matter particles must be (almost) stable on cosmological timescales. In many models 
this stability is provided by a conserved Z 2 quantum number generically called parity. In 
supersymmetric models its role is played by the i?-parity, while in more phenomenologically 
oriented models such as the Inert doublet model, j2], or the minimal singlet model, j3], the 
Z 2 symmetry is imposed by hand when constructing the lagrangian. In these models, all the 
SM particles including the SM-like Higgs boson have positive parity, while the dark sector 
particles are of negative parity. The lightest among these negative parity particles, which we 
will generically denote as d, is stable and represents the dark matter candidate. 

It is natural to ask whether dark matter can be stabilized by a conserving discrete quantum 
number taking values in a group other than Z 2 . This idea was explored in a descent number of 
works, in which both abelian [U [5j [6] and non-abelian [7] finite groups were usecu Requiring 
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1 If one takes seriously the argument that quantum- gravitational effects violate any global discrete symmetry, 



that d transforms non-trivially under the discrete group implies that it cannot be a truly 
neutral particle. It also prohibits direct two-particle annihilation dd — y X S m, where X S m 
is any set of SM particles. However more complicated processes involving several ef s can 
take place. This certainly changes the kinetics of dark matter evolution in the early Universe 
and its relic abundance after the freeze-out. If in addition one assumes that an asymmetry 
between d and its antiparticle d* is generated at a high energy scale in a way similar to the 
baryon asymmetry of the usual matter [9] , then the present day behavior of this asymmetric 
dark matter is not dominated anymore by dd* annihilation and can lead to characteristic 
observational signatures. 

One particular class of groups used to stabilize dark matter are cyclic groups Z p , see 
jU |5] and references therein. With this choice, all fields are characterized by a conserved 
quantum number which we will call the Z p charge q and which is additive modulo p. The 
usual assignment is that all the SM fields including the SM-like Higgs boson have q = 0, while 
the dark matter candidates have nonzero Z p charges. This opens up the possibility of novel 
two-particle processes such as dd — y d*XsM, which was called semi-annihilation in [TO], or 
even multiparticle versions of it, "multiple semi-annihilation". Besides, if the model allows 
for existence of several dark matter candidates di with different charges <&, then inelastic two- 
particle processes didj — y d k X SM are also possible. Such processes, too, have impact on the 
kinetics of the dark mater abundances in the early Universe, [11] . 

Examples of Z p -stabilized dark sectors often involve a variety of new fields which interact 
via lagrangians designed specifically to incorporate a given symmetry group, see e.g. [5]. Even 
when linked to electroweak physics, these models involve new fields with different electroweak 
quantum numbers (doublets, singlets, etc.), [IT] . Indeed, if one assumes that extra fields come 
from a hidden sector coupled to the SM fields via "portal" operators [12], then they must be 
electroweak singlets by construction. 

In this paper we demonstrate that Z p -stabilized scalar dark matter can easily arise in multi- 
Higgs-doublet models. This, perhaps, is not surprising on its own. A less trivial fact is that 
even with few doublets one can get Z p with a rather large p. To be precise, in models with N 
Higgs doublets, Z p with any p < 2 Ar_1 is realizable in the scalar sector. We will show that this 
fact can be instrumental in avoiding semi- annihilation processes even with multi-component 
dark matter sectors. 

The structure of the paper is the following. In Section |2] we quickly review what is known 
about symmetries in the scalar sector of iV-Higgs-doublet models. Then, in Section [3] we 
give a 3HDM example with dark matter candidates stabilized by the Z3 symmetry group. 
In Section H] we show what is possible with four doublets and consider in some detail the 
Z7-symmetric 4HDM, in which one can avoid semi-annihilation processes. We end the paper 
with a discussion and conclusions. 

2 Scalar sector of NHDM and its symmetries 

The iV-Higgs-doublet model is a conceptually simple extension of the SM Higgs mechanism. 
It is driven by the idea that the Higgs fields, similarly to fermions, can come in several 
generations. Its simplest version with only two doublets, 2HDM, is of special interest because 



one should require that this quantum number arises as a remnant of a C/(l) gauge symmetry spontaneously 
broken at a high energy scale, [5]. 



it mimics the Higgs sector of MSSM, and it has been thoroughly studied in the last four 
decades, see [13] and references therein. In addition, many particular models employing more 
than two scalar doublets have also been proposed [H]. 

In NHDM one introduces N doublets of complex scalar fields <f>{, i = 1, . . . , N, with elec- 
troweak isospin Y = 1/2, and constructs the self-interaction renormalizable Higgs potential 

V = YMh) + Z ljkl (4<f>jM<f>i) , (1) 

where all indices run from 1 to N. The free parameters of the potential are written as com- 
ponents of tensors Y^ and Z^i] in the most general case there are N 2 (N 2 + 3)/2 independent 
parameters. Once these coefficients and the Yukawa couplings are provided, the model is com- 
pletely defined and the entire phenomenology should follow. In practice, however, inferring 
these consequences directly from the lagrangian is impeded by algebraic obstacles at the very 
first step, namely, the minimization of a sufficiently generic potential. The consequence is that 
only very few general results are known for N > 2, [T5l ITBl [T7] . 

One particular issue which is of much importance and where certain progress has been 
recently made concerns accidental symmetries which can be encoded in the scalar sector of 
NHDM. These are transformations that mix several Higgs doublets (also called Higgs-basis 
transformations) but still leave the potential invariant due to specific patterns in the tensors 
Yij and Z^i- Although particular models with several Higgs doublets based on various sym- 
metry groups have been proposed and studied over the last decades, no attempt at systematic 
classification of possible symmetries was made until very recently. 

One of the problems here comes from the observation that just imposing a certain symme- 
try group G on the Higgs potential can often lead to potentials symmetric under a larger group, 
which includes G as a subgroup. This feature was discussed in [16] where one-parametric sym- 
metry groups were studied for 3HDM and in [17] where an attempt to understand symmetries 
in NHDM via geometric constructions in the space of bilinears was made. Therefore when 
describing symmetries of the model, one should focus only on the true symmetry groups, the 
ones which are not automatically extended to larger groups. These groups were called in [17] 
"realizable groups". 

Unfortunately, no systematic way to reconstruct the true symmetry group of a given NHDM 
potential for iV > 2 is known so far. One even does not know the list of realizable symme- 
try groups possible for N = 3. However a step forward was recently made in [18], which 
completely characterized all abelian groups of Higgs-family transformations and generalized 
CP-transformations that can be realized as symmetry groups of the scalar sector of NHDM. 
In particular, in what concerns realizable cyclic groups, it was proved there that for a model 
with N doublets one can construct the Higgs potential symmetric under the group Z p with any 
1 < p < 2 W_1 (for N = 3 this conclusion was known before, [16]). One can therefore wonder 
if these models can be used to construct dark matter sectors stabilized by these symmetry 
groups. 

An electroweak symmetry breaking (EWSB) model with Z p -stabilized scalar dark matter 
must satisfy several conditions. First, the entire lagrangian and not only the Higgs potential 
must be Z p -symmetric. The simplest way to achieve this is to set the Z p charges of all the SM 
particles to zero and to require that only one Higgs doublet (the SM-like doublet) couples to 
fermions. The Z p charge of this doublet must be zero, and it does not matter which doublet 
is chosen to be SM-like due to the freedom to simultaneously shift the Z p charges of all the 
doublets. Second, the Z p symmetry must remain after EWSB. This is possible when only the 



SM-like doublet acquires a non-zero vacuum expectation value (v.e.v.). Third, if we insist on 
Zp-stabilization, that is, we require that not only decays but also 2-, 3-, . . . , (p — l)-particle 
annihilation to SM fields are forbidden by quantum numbers, then the dark matter candidates 
must have Z p charge q which is coprime with p. 

Below we show that multi-Higgs-doublet models can easily satisfy these conditions. We 
start first with the simplest model of this kind, Z3-symmetric 3HDM, and then we outline the 
plethora of kinds of microscopic dark matter dynamics possible in models with four doublets. 

3 Z 3 -symmetric 3HDM 

A natural way to implement the Z3 symmetry in 3HDM would be to construct a potential 
symmetric under 0i — > 2 — >• 03 - *■ 0i- However, upon an appropriate Higgs basis change, this 
transformation will turn into pure phase rotations of certain doublets. In fact, it can be proved 
that any abelian subgroup of SU(N) can be mapped onto a group of (possibly correlated) phase 
rotations of individual doublets, see [18] for explicit construction. Therefore, we will always 
use below the phase rotation representations of the cyclic symmetry groups. Also, to keep the 
notation short, we will describe any such transformation 0, — > e iai <fii by providing the iV-tuple 
of phases a«. 

A scalar potential invariant under a certain group G of phase rotations can be written as 
a sum V = Vq + Vg, where V is invariant under any phase rotation, while Vq is a collection of 
extra terms which realize the chosen symmetry group. The generic phase rotation invariant 
part has form 

V o = E [~ m M<t>i) + A M (0l0,) 2 ] + J2 [\0Mi4>]^) + KMh)^)} > ( 2 ) 

i ij 

while Vq obviously depends on the group. In particular, for the group Z 3 in 3HDM we have 
Vi 3 = A 1 (0 3 1 )(40 1 ) + A 2 (0i0 2 )(0 t 3 2 ) + A 3 (0 2 3 )(0i03) + h.c, (3) 

where at least two of the coefficients Ai, A 2 , A 3 are non-zero (otherwise, the potential would 
have a continuous symmetry). This potential is symmetric under the phase rotations generated 
by 

« = y(0, 1, 2), a 3 = l, (4) 

In fact, the assignment of these charges to the three doublets is completely arbitrary, and 
the group generated by the generator a with permuted charges has the same action of the 
potential. In (j4]) we simply chose 0i to be the SM-like doublet. 

Whether this symmetry is conserved or spontaneously broken depends on the pattern of 
the vacuum expectation values. If we insist on conservation of the Z3 symmetry, we require 
that (4>1) = V\ = v/y/2, (02 1 ) — (0°) — 0. There is nothing suprising that the potential 
Vq + Vz 3 can have a Z 3 -symmetric global minimum upon an appropriate choice of coefficients. 
The question is whether it requires any fine-tuning or not. Below we show that it does not, 
and the minimum of the type (v 1, 0, 0) arises in a sizable part of the entire available parameter 
space of this model. 

First, we note that if (v~i, 0, 0) is an extremum of Vq, then it is also an extremum of 
Vq + Vz 3 because the extra terms contain X only linearly and quadratically. Therefore, when 



constructing a model example, one can first build Vq with a global minimum at (t>i, 0, 0) and 

then add a sufficiently weak Vz 3 so that this point remains a minimum. 

Now, turning to minimization of Vo, suppose that we search for the neutral minimum with 

a generic complex v.e.v. pattern (vi, v 2 , ^3). By introducing pi = \vi\ 2 > 0, we can rewrite Vo 

as 

1 1 

V = -Mip-i + -AijPiPj = -Aijipi - Pi){pj - Pj) + const . (5) 

Here Mi = (ml, m 2 , m§), Ay is constructed from Ay and A„- in an obvious way, and pi = 
(A i j)~ 1 Mj. Positivity condition on the potential guarantees that Ay is a positive definite 
matrix, therefore its inverse exists (we omit here the degenerate situations when the quartic 
potential has flat directions). 

The allowed values of the pi populate the first octant (pi > 0) in the three-dimensional 
euclidean space. Due to (jSJ), the search for the global minimum can be reformulated as the 
search for the point in the first octant which lies closest to pi in the euclidean metric defined 
by Ay. Clearly, if pi itself lies inside the first octant (p\, p 2 , P3 > 0), then the global minimum 
is at pi = pi. If pi lies outside the first octant, then the closest point lies either on a face, or on 
an edge of the first octant, or at the origin. The global minimum is unique by the convexity 
arguments. It is now clear from this geometric construction that the entire space of all possible 
vectors pi can be broken into several regions of non-zero measure which correspond to all of 
these possibilities. In particular, the region corresponding to vacuum alignment of the type 
(t>i, 0, 0) also has a non-zero measure and fills a sizable part of the parameter space. In this 
sense, this vacuum pattern does not require any fine-tuning of the coefficients of the potential. 

Having established that the required vacuum pattern is generically possible, we now switch 
to a simple version of the model. This is done only to simplify the presentation of the argument; 
if needed, the calculations can be repeated for a generic potential. Namely, we now choose 
mf > 0, 7713, m 3 < an d take Ay = 2X 5ij, so that the potential becomes 

V = -m?(0t0i) + \rr%\(4<h) + \m 2 3 \(4fo) + ^0 [{<f>\<t>if + {^Ifa? + {<&fo? 

+A 1 (0 3 1 )(0 2 1 ) + A 2 (0l0 2 )(0^ 2 ) + A 3 (0 2 03)(0l03) + h.c. (6) 

By construction, its global minimum is at (0°) = (v/y/2, 0, 0), where v 2 = m1/\ . In order 
to find the mass matrices, we write the doublets as 
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G+ 
±(v + h + iG°)j ■ yz 2 j ■ yz 3 



(7) 



Here h is the SM-like Higgs boson, G° and G + are the would-be Goldstone bosons, while 
w 2 , W3 and z 2 , z 3 are charged and neutral Higgs bosons, respectively. Fields w^ and 2:2,3 
have well-defined Z 3 -charges: q = 1 for w^, z 2 as well as for w 3 , z^ and q = 2 ( = — 1 mod 
3) for w 3 , z 3 and w 2 , z 2 . Note that in contrast to the usual practice, we describe the neutral 
Higgs bosons in the second and third doublet by complex fields rather than pairs of neutral 
fields. The reason is that, by construction, the fields corresponding to the real and imaginary 
pairs of z's have identical masses and coupling constants. In any process that can arise in this 
model, these two fields are emitted and absorbed simultaneously, so they can be described by 
a single complex field Zi. 

The SM-like Higgs boson has mass mj* = 2ml, while the masses of the charged Higgs 
bosons are m 2 ± = \m 2 \, m 2 ± = \m 2 \. Neutrals with equal q can mix, which indeed happens 
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at Ai 7^ 0. We descibe the resulting mass eigenstates by complex fields d and D {m d < m D ) 
both having q — 1: 

d = cos a z 2 + sin ae~ ! ^ z% , £> = — sin ae l/3 z 2 + cos a Z3 , 

tan 2a = — — -. — T — -. — 5-7, /3 = argAi 
A I ml 



m 




^- — 2 — ^yH-w + yX- (8) 



N 2 



Note that within this model we have 

m d < "V , rn w ± < m D , m 2 ± + m 2 ± =m 2 d + m 2 D . (9) 

Z o 2 3 

The triple and quartic interaction terms arising from Vq and Vz 3 specify the dynamics of the 
dark matter candidates. The lightest particle from the second and third Higgs generations is d, 
and it is stabilized against decaying into the SM particles by the Z 3 symmetry. As for heavier 
particles, triple interactions lead to their decays such as D — > dh, D — > dZ, and D — > w^W - , 
D —7- w^W + if allowed kinematically. If the mass splitting between the d and D is small, 
then these processes involve virtual h, Z, etc. which then decay into the SM particles. In this 
aspect, D decays are similar to weak decays. Charged Higgs bosons w 23 will also decay to d 
or d* plus SM particles. 

In the case of symmetric dark matter, the main process leading to depletion of dark matter 
after electroweak symmetry breaking is the direct annihilation dd* —¥ Xsm, and the semi- 
annihilation reaction discussed below is only a correction to this process. Still, it might be 
possible that this correction leads to a sizable departure of the kinetics of the dark matter in 
the early Universe and affects the relic abundance at the freeze-out, [TTj . 

The situation becomes more interesting in models of asymmetric dark matter, [9], in which 
an asymmetry between d and d* is generated at a higher energy scale. It is possible for example 
that upon electroweak phase transition almost all d's annihilate with d* into the SM sector, 
leaving behind a certain concentration of dark matter candidates d. In the present epoch, d 
can scatter elastically, dd —¥ dd with a e i, oc A 2 ,, but they can also initiate semi-annihilation 
processes such as dd — > d*X SM with a subsequent annihilation of d* with a d. This possibility 
originates from the following quartic terms in the scalar potential 

—=hddd cos a sin ae ll3 {\2 cos a + A3 sin ae 1 ^) + h.c. (10) 

v2 

Depending on A's, this process can be as efficiently as the direct annihilation in the usual 
annihilating dark matter models, or it can be suppressed by the small coupling constant. 

Finally, the same interaction terms also generate the triple annihilation processes ddd — > 
h — > Xsm, whose rate is, however, suppressed at small densities with respect to the semi- 
annihilation. 

4 Avoiding semi-annihilation 

The presence of the hddd terms in the interaction lagrangian in the previous example, which 
were responsible for the two-particle semi-annihilation process, was due to the Z3 symmetry 
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group. One can wonder whether two-particle semi- annihilation can be avoided by employing 
a Z p group with larger p. In this present section we show that it is indeed possible in a model 
with four Higgs doublets. 

According to [18], one can encode in the 4HDM scalar sector any group Z p with p < 8. 
Note that in order to avoid continuous symmetry, one must accompany the phase-symmetric 
part of the potential Vo with at least three distinct terms transforming non-trivially under 
phase rotations. In Table 1 we give the list of these symmetry groups together with examples 
of the three interaction terms and the phase rotations that generate the corresponding group. 



group interaction terms 



~% (it2), (its), (iuy 

Z 3 (3+2), (It3)(4t3), (lU)(lt2) 

Z 4 (3+2), (1+3X4+3), (1+4) 2 

Z 5 (4+3)(2+3), (3+2)(l+2), (4+l)(3+l) 

Z 6 (4+3)(2+3), (3+2)(l+2), (1+4) 2 

Z 7 (4+l)(3+l), (4+3)(2+3), (4+2)(l+2) 

Z 8 (4+3)(2+3), (4+2)(l+2), (1+4) 2 



phase rotations 


f (0, o, 


0. 


1) 


f (0, 1, 


1. 


2) 


f (o, i, 


1- 


2) 


T(°. X ' 


2. 


3) 


f (o, i, 


2. 


3) 


f (0, 2, 


3. 


4) 


f (0, 2, 


3, 


4) 



Table 1: Cyclic groups realizable as symmetry groups in the scalar sector of 4HDM; (a+fe) is a 
short notation for (</>+</>&). 



Two remarks concerning this table are in order. First, we stress that all these groups are 
realizable, so that if the three terms in each line are written down with non-zero coefficients, 
there is no phase transformation other than multiple of the generator and the common overall 
phase shift that leaves them invariant. Second, since the potential is invariant under the 
common phase shift of all doublets, one can freely add additional equal phases to the generators 
shown in the third column and, in addition, one can permute the doublets. For example, the 
last line of this table can be replaced by (3+2)(l+2), (3+l)(4+l), (4+3) 2 , which is symmetric 
under the Z 8 group generated by phase rotations ^ (0, 1, 2, 6). 

The patterns of phase shifts given in this table allow for construction of various dark sectors 
with different possibilities for dark matter dynamics. Here, we do not aim at a complete 
classification of these possibilities but would like only to show that there are examples in 
4HDM in which semi-annihilation of dark matter candidates is also forbidden. 

To this end, let us consider the Zy-symmetric 4HDM with the potential V = Vq + Vz 7 with 

v Z7 = \i(<A<h.)(4<l>i) + HtplfoMfo) + H4>lfoK4fo) + h.c. (ii) 

As before, we assume that only the first doublet couples to fermions. Using the technique 
of the previous section, one can easily construct the potential Vq with the global minimum 
at {v/y/2, 0, 0, 0). Then expanding the doublets similarly to ([7]), one observes that z 3 and 
z\ have Z 7 -charges q = 3 and mix via the Ai term leading to mass eigenstates d and D. In 
addition, we have the field Z2 with Zy-charge q = 2 and electrically charged Higgs bosons wf 3 4 
with appropriate Z 7 -charges. 



By adjusting free parameters, one can easily make d the lightest among the particles that 
transform non-trivially under Z 7 . Then, the other particles will either eventually decay to d or 
d* plus SM particles or will be stable representing an additional contribution to dark matter. 
Again, if the asymmetry between d or d* exists and if the rate of their annihilation is high, 
then after the freeze-out we are left with the gas predominantly made of d's. 

The subsequent microscopic dynamics depends on the interactions between d's and Z2's 
which follow from (ITT|) . The relevant terms are dz^z^ times a SM field from the A2 term and 
dddz\ from the A 3 term. One- or two-particle processes such as d — > z^z^, dd — > d*z 2 , and 
dd — > Z2Z2Z2 are all kinematically forbidden. Multiple collision kinetics depends on whether 
m Z2 < 3md or not. If Z2 is not too heavy, then the "triple semi-annihilation", ddd —¥ z%Xsm, is 
kinematically allowed and will create a population of Z2 even if it were absent before. However, 
z 2 will get depleted by the semi-annihilation process Z2Z2 — > d*XsM- So, if one starts with a 
certain concentration of Z2, d and their antiparticles, then 22 will die off at a higher rate than 
d's. In stationary conditions, the terminal concentrations will be those equilibrating the rates 
of the following Qd tree- level scattering with intermediate 22 's: 

6d^z 2 z 2 XsM^d*X' SM , (12) 

and the subsequent annihilation of d* . The net result of this chain will be the "7c?-burning 
process", Id — > Xsm, the bottleneck in this chain being the triple-d process ddd — > z%Xsm- 

On the other hand, if m Z2 > 3rrid, then ddd — > Z2 is kinematically forbidden, while the 
inverse process leads to a quick Z2 decay. In this case, one can still burn d's via the tree-level 
process with intermediate virtual Z2 's: 

dddd —¥ d*d*Z2Z2 — > d*d*d*XsM 1 (13) 

The net result will be the same 7(i-burning, but the bottleneck process is now the Ad collision, 
whose rate is even stronger suppressed. 

5 Discussion and conclusions 

The main purpose of this paper is to demonstrate that there already exists a phenomenological 
template for scalar dark matter models stabilized by cyclic groups larger than Z 2 . This 
template uses several electroweak Higgs doublets decoupled from fermions, and it represents 
one of the simplest extensions of the Standard Model. Remarkably, models with few doublets 
can easily accommodate dark sectors which are stabilized by a large list of discrete groups and 
which display various kinds of microscopic dynamics. In particular, we gave explicit examples 
of dark sectors where the bottleneck process leading to depletion of asymmetric dark matter 
can be a 2-particle, 3-particle or 4-particle semi-annihilation. We stress that these models do 
not require any serious fine-tuning. We only ask for the presence of terms invariant under the 
chosen symmetry group but do not constrain coefficients in front of these terms. 

In certain aspects these models resemble the Inert Doublet Model, [2J, but in the other 
they rely on symmetry patterns that arise only with several doublets. In this respect, such 
models can be viewed as "multi-inert" doublet models although this name of course does not 
completely specify the microscopic dynamics. 

Exploring the observational consequences of each sort of microscopic dynamics is a separate 
task. It should include study of the dark matter kinetics in two situations. First, one obviously 
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needs to track down the dark matter evolution in the expanding Unverse after the electroweak 
phase transition and determine the freeze-out abundances. Analysis of [H] already proves that 
semi-annihilation processes can be important, but it remains to be understood how sensitive 
the evolution is to the exact microsopic dynamics. 

We would like to stress that studying this problem in the context of multi-Higgs-doublet 
models can be much subtler than it looks at first glance due to multiple phase transitions 
near and below the electroweak scale. Indeed, even in the two-Higgs-doublet model a single 
electroweak phase transition can split into a sequence of several phase transitions of different 
nature, both in the general case [19] and in the Inert doublet model [20]. One can expect that 
even longer chains of phase transitions can be possible in multi-doublet models. Note that the 
last among these phase transition can in principle happen at temperatures much lower than 
the nominal electroweak temperature scale. Consequently, the Universe might have evolved 
through a sequence of vacua with different, and perhaps exotic, properties. Phase transitions 
between these phases could have led to complete restructuring of the particle mass spectrum, 
both within the SM and in the dark sector; particles which are stable in one phase can be 
unstable in another. All these delicate details, as well as the thermodynamics of the phase 
transitions themselves, can modify the evolution of the dark sector. None of the existing 
evolution codes can adequately address these intricacies. 

The second situation where the microscopic dark matter dynamics can make an impact 
is the present epoch evolution at astrophysical sites of elevated dark matter concentrations 
(galactic centers, interiors of compact stars, etc., [21] ). Since multi-particle processes are 
involved, the sensitivity to the dark matter density will be different from that of the usual 
two-particle annihilating or with truly non-annihilating dark matter. 

For example, it is known that dark matter with sufficient elastic cross section can get 
captured inside neutron stars, [21]. In models with asymmetric truly non-annihilating scalar 
dark matter, its accumulation can lead to formation of the Bose-Einstein condensate (which 
means that the occupation number in the phase space can become large) or even to collapse in 
a tiny black hole [22]. In a certain region of parameter space, this black hole will destroy the 
host neutron star sufficiently quickly compared to the typical neutron star lifetime; therefore, 
this region is excluded by observations. In the case of multi-inert dark matter accumulated 
inside a neutron star, multiple annihilation processes will effectively enter the game as the 
density reaches a certain threshold, precluding black hole formation and avoiding the above 
constraints. 

In conclusion, we showed that multi-Higgs-doublet models can naturally accommodate 
scalar dark matter candidates protected by the group Z p . For a model with iV doublets, 
the values of p can be as large as 2 N ~ l . These models do not require any significant fine- 
tuning and can lead to a variety of forms of microscopic dynamics among the dark matter 
candidates (allowing or forbidding semi-annihilation, offering different routes to multi-particle 
annihilation, etc.). 
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